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1 Introduction 

In this paper we study existence and uniqueness of bounded solutions to 
nonlocal nonlinear initial- value problems of the following type: 

P d ^ u +(-^y V] = o xeR,t>o (L1) 

u = uq x G K, t = 0; 

here p = p(x), usually referred to as a variable density, is a positive func- 
tion only depending on the spatial variable x, {^—-j^^j denotes the one- 
dimensional fractional Laplacian of order 1/2, m > 1 and uq is a nonnegative 
and bounded function. 

The interest in the study of nonlocal problems has grown significatively 
in the last years and the analysis of nonlinear equations involving fractional 
powers of the Laplacian has become an area of intense research. Such prob- 
lems in fact arise in many physical situations (see, e.g. pQ, |18j , |19j ) in 
particular in the analysis of long-range or anomalous diffusions. From a prob- 
abilistic point of view, the fractional Laplacian is the infinitesimal generator 
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of a Levy process (see [2]). Then, losely speaking, uniqueness for problem 
with m = 1 corresponds to the fact that the Levy process associated 

to the operator j -§^x\ > starting from any point in K, does not attain 
infinity. 

If the nonlocal operator (— in (jl.ip is replaced by the classical 

second order partial derivatives — -§-^1 we get 



■^- = ,el,oo (L2) 
i e 1, ( = 0, 

which is the initial value problem for the porous medium equation with vari- 
able density in one space dimension. Problem (|1.2p and its counterpart in 
R N , that is, 

(pd t u + Au m = xeR N ,t>0 
\u = uq x eR N ,t = Q, 




occurs in various situations of physical interest (see, e.g., [53]) and have been 
largely investigated in the literature during the last two decades (see [14], 
[15] . [Id] . [17] . [2D] , [2T] , [24] , [28] -[25]). More recently, similar problems on 
Ricmannian manifolds have been addressed in [25]-[27j. 

It is well-known that the asymptotic behavior of the varying density p may 
influence uniqueness of solutions for problem (|1.3[) . Briefly, for N — 1 or N = 
2 problem (|1.3[) is well-posed in the class of bounded solutions not satisfying 
any additional conditions at infinity (see [T7]), provided p G L co (R N ). The 
situation is different if N > 3. Indeed, in this case uniqueness prevails when 
the density p goes to zero slowly as |x| — > oo. If instead p goes to zero fast as 
| a; | diverges, then nonuniquencss of bounded solutions occurs (see [T5], [TP] , 

Concerning the regularity assumed for initial data and solutions we point 
out that in [20] . |28j-|29| weak energy solutions are dealt with, and uq is 
supposed to belong to Lp(R w ), the space of measurable nonnegative functions 
/ satisfying J RN f pdx < oo. In [T5], [TDJ [2T] and [23], instead, bounded 
initial data and so-called very weak solutions are treated. One of the main 
issues in dealing with very weak solutions lies in the fact that, in general, 
these solutions do not have finite energy in the whole R*. 

Very recently, in [11| , the nonlocal nonlinear initial- value problem 

f«9 t u + (-A)f [u m ] =0 x€R N ,t>0 
lit = u Q x g M w , t = 



(0 < a < 2) has been studied. The particular case a = 1 has been considered 
in jlOj ; more precisely, existence, uniqueness and properties of weak solutions 
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to (|1.4p have been established assuming uq G L 1 (R Ar ). Furthermore, problem 

d t u+ (-^Y pog(l + u)] =0 i£l,i>0 
u = uo x G K, t = 

has been addressed in [12] . 

The motivation in studying problem (jl.ljl is twofold. In fact, on one hand, 
when p = 1, (|1.1[) becomes problem (|1.4j) with a = N ~ 1, thus problem 
can be regarded as a generalization of problem (|1 .4[) to the case in which a 
variable density p is taken into account. On the other hand, problem 
can also be considered as nonlocal version of problem (|1.2[) . 



The aim of this paper is to prove well-poscdncss of problem in the 
class of bounded solutions not satisfying any additional conditions at in- 
finity. We shall consider bounded initial data uq and, consequently, very 
weak bounded solutions (see Definition 12. 1|) . Let us mention that our re- 
sults differ from those in |10| where p is assumed to be identically equal to 
1. Moreover, also in the case p = 1 considered in |10| uniqueness of very 
weak solutions is not proved, while uniqueness of weak energy solutions is 
shown (see [101 Lemma 4.1]). Very weak solutions to problem (| 1 .4[) in the 
space C([0, oo); L^(R N )) (ip being a suitable weight decaying at infinity) 
have been considered in [3] where uniqueness of such type of solutions is also 
shown provided < m < 1. 

To the best of our knowledge, the uniqueness result presented in this 
paper is new also in the linear case m = 1, that is for problem 

pdtU+ {-&Y M= ° x&R > t>0 (1.5) 

U = Uq X G K, t = 0. 

Some results for nonlocal linear parabolic equation with a variable density 
have been established in [7], but not for problem (|1.5[) . 

One can expect that if N > 2 the nonlocal problem 

{pd t u+(-A)i [u m ] =0 x G M. N , t > 
\u = u xeR N ,t = 

may exhibit uniqueness or nonuniqucness of solutions depending on the be- 
havior at infinity of the density p. This would reflect in the nonlocal case 
the the situation described above for problem (|1.3[) . However, some tech- 
nical difficulties prevent us from adapting the methods used in this paper 
to investigate uniqueness of very weak solutions in any space dimension; see 
Remark 14.41 We postpone to a forthcoming paper the proof of existence and 
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uniqueness of solutions to problem 

ipd t u+ (-A)f [u m ] = xeR N ,t>0 
\u = u x£R N ,t = 

for any < a < 2 and any N > 1. In this case, according with [TD]-[TT] , we 
shall consider weak energy solutions and completely different arguments. 

The paper is organized as follows. In Section [2] we introduce the setting 
of the problem and give the precise definition of solutions we are dealing 
with. We also recall some basic facts about fractional Laplacian, mainly its 
realization trough the harmonic extension, which permits to look at problem 
p. II) in a "local way" as a quasi-stationary problem with dynamical boundary 
conditions. In Section [3] we prove existence of solutions to problem see 
Theorem 13.11 Even if the general strategy of the proofs in this Section goes 
along the same lines as in [TO]- [TTj , there arc some differences here that will 
be expedient as we will point out in Remark 13.81 We should remark that all 
results in Section [3] are valid in general for any N > 1. In subsequent Section 
U]we establish uniqueness of solutions; see Theorem 14.11 The proof exploits 
the above mentioned realization of the fractional Laplacian combined with 
properties of a suitable family of test functions we construct in Lemma 14.21 
Here the restriction N = 1 will be crucial as explained in Remarks 14.41 and 
1431 



To conclude, let us mention that our techniques also applies to more 
general fractional nonlinear diffusion equations with variable density, such as 

pd t u+(--^y [G{u)} = xeR,t>0, 

provided G satisfies suitable conditions as in [15], [24]. Moreover, changing 
sign bounded solutions could also be considered. However, for sake of sim- 
plicity, we limit ourselves to the case of nonnegative solutions and G(u) = u m 
(to > 1). 



2 Problem setting and assumptions 

i 

We recall that the nonlocal operator ^— is defined for any function (p 

belonging to the Schwartz class by 

*)%_Ip.v./2feL^U. 

ox ) I" Jr \x-y\ 

1 

As is well-known, (^-g^j can be also defined in many other ways; we refer 
the reader to [1], [B], [T3], for a comprehensive account on the subject. In 
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particular, if tp is a smooth and bounded function defined in 
consider its harmonic extension v to the upper half-space 



fi := Rj. 

Hence (see g], [BJ, [TO], [12]), 

9u(a;, 0) 
<9y 



{(a:, j/): x G R, y > 0}. 



9x 2 



for all x G 



In the sequel, we always make the following assumption: 
f(i) p G C(R), /? > in R; 

1 (ii) uq is nonnegative, continuous and bounded in 



(Ac) 



Following [TO] and [TT], using the harmonic extension, we shall rewrite 
the problem (jl.ip in terms of local differential operators. More precisely, 
solving problem (jl.ip is equivalent to solve the quasi-stationary problem with 
a dynamical boundary conditions: 



Aw = 

dw 

oy 



at 



(x,y) G fi, t > 

i g r, t > o 

i 6 f, t = 0; 



(2.1) 



here r := n = 0} = R. 

Now, take any open bounded subset D C with n smooth and 
3D FIT ^ 0. Then the exterior normal to 3D exists almost everywhere; we 
denote such vector by v. Let T > 0,ip G C 2 x ' l t {D x [0,T]), ^ > 0, ip = 
in (5D\r) x [0,T], V ^ in (9L> n T) x [0,T]. Formally, multiplying the 
differential equation in (|2. 1[) by ^, integrating by parts and considering the 
initial and the dynamical boundary condition we get: 



= 



J D 



tp Aw dx dy dt 



(V^Vw) dx dy dt 



J D 
T 



J D 



(Vip, Vw) dx dy dt - 



, 9w 

W ttz dx dy dt 
dv 

i ® w i i 
tp -jr— dx dt 

o JdDnr oy 



o J D 



w Aip dx dy dt — 



T 



w-^z dS dt 

o JdDnn ® v 



w— — dx dt 

o JdDnr vy 



pipdtudxdt 



o JdDnr 
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/ / w Aip dx dy dt - / / 
Jo Jd Jo JdDnn 



w-^t dS dt 
av 



f T f dtp , , 
/ / w— dim 

Jo JdDnr vy 



pudtipdx dt 



IdDnr 



lo JdDnr 
p(x) u(x, T)ip(x, 0, T) — uo(x)ip(x, 0, 0) 



dx. 



Given a function / e W,'(fl) we denote by /|r its trace on T, which is 
in Lf oc (T). In view of previous equalities we give next 

Definition 2.1. A solution to problem (|2.I[) is a pair of functions (u,w) 
such that 

• it 6 L°°(T x (0,oo)), u > 0; 

• w <E L°°(fl x (0,oo)), i«>0, Vwe Lf oc (Ti x (0,oo)); 

• w|rx(o,oo) = u m ; 

• for any T > 0, D and ip as in (??) there holds 



T 



o Jd 



w Aip dx dy dt 



o JdDnn 



w-—z dS at 

av 



+ 



u^dxdt 
o JdDnr °V 



T 



pudtipdx dt 



Jo JdDnr 

-J p(x) u{x,T)ip{x 7 0,T) - u (x)(p(x,0,0) dx = 0. (2.2) 

Here and hereafter we say that / £ Lf oc (fl x (0, oo) when / G L 2 (fl' x 
(0, oo)) for every open subset fl' CC ft. Note that, since Vw E Lf oc (f2 x 
(0, oo)), the equality u>|rx(o,oo) = u m makes sense. 



3 Existence of solutions 

In this Section we prove existence and some properties of solutions to problem 

CD)- 

Theorem 3.1. Let assumption flApD &e satisfied. Then there exists a solution 
(u,w) to problem (|2.ip . 

The proof of Theorem 13.11 is divided in several steps. Let us fix some 
notations. For any R > we denote by B R = B N (0,R) the Euclidean open 
ball in 1*; we set 

fin := fin B 2 R , 

T R := dfl R n {y = 0} = , 

S fl := n {y > 0}. 
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If / £ W /1 ' 2 (r2fl), we denote by f\r R ,f\s R the trace of / on and S_r, 
respectively. 

For any fixed e > 0, R > and g G L°°(Tr) we consider the following 
auxiliary elliptic problem in D,r: 

'Aw fl = (x,y)en R 

+ p(vr)™ = pg lEffi (3.1) 
= (x,y) G 

Definition 3.2. A solution to problem (|3.1[) is a pair of functions (zr,vr) 
such that 

• z R 6 i°°(r fl ); 

• for any tp G C (CIr), ip = on E/j i/iere ftoWs 



/ (Vufl, V<p)dxdy + - f p(z R -g)<pdx 
Jn R £ Jr R 



(3.2) 



In the following Lemma we establish existence, uniqueness and some prop- 
erties of solutions to problem (|3.ip . 

Lemma 3.3. Let R > and e > be fixed. For any g G L°°(T), there exists 
a unique solution (zr,vr) to the problem ()3.1j) . Furthermore the following 
properties hold: 



i. Contr activity: the mapping 



1 

9 i-> -z R 
P 



is a contraction in the norm of L}(Yr), namely, if vr and Vr are 
solutions of (|3.1[) corresponding to g and g respectively, then 

p{z R - z R ) + dx < / p[g-g], dx; (3.3) 



in particular, if g > then zr>0 in £Ir; 
ii. Uniform boundedness: for any e > and any R > 

< vr < ||<?||£o (r) m n R , (3.4) 

provided g > 0; 
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Hi. Monotonicity: if R' > R then vr/ > vr in Qr, provided g > 0. 

Proof. Let H 1 ^) := {/ G H 1 ^) : f\s R = 0}. Existence of solutions 
easily follows by standard arguments, since the functional J : H 1 (Hr) — > M 
defined by 



1 



J( v ) := o / \Vv\ 2 dxdy-i / pv ™ dx — I pvgdx 



2 Jn R m + 1 J Tn 



is coercive in in fact, for suitable constants G\ = C\{R) > and 

C 2 = C 2 {R) > 0, we have 

J(v) > — f |Vw| 2 da;dw+ ( minp| / v m dx 

yJ -2Jcin V r « Jm + lJ rR 

- y^pj J r v 9 dx - G i\\ v \\m(a R ) - C 2 \\v\\ H i {nR) . 

Now, to show contractivity, take a smooth monotone approximation p of 
the sign function, and set ip(x) := p[vr(x) — vr(x)], x G FIr. Clearly, such 
a tp can be used as a test function in the weak formulation of Definition 13.21 
We get 

= / p'{v R - vr)\V(v r - v R )\ 2 dxdy 



Jn R 

+ - [ p( ZR -~z R )p[{z R ) m -{z R ) m ]dx-- [ p(g~~g)p[{zR) m -(~ZR) m ]dx. 

e JTn e JTr 

Passing to the limit in the approximation p, we discover 

p(z R - zr) + dx < / p {g - g) sgn(z]£ - z R ) dx 
r R Jr R 



which immediately gives (|3 . 3|) . Furthermore, if g > 0, (J3T3J) implies zr > 
in Tr. From (|3.3[) uniqueness can be easily deduced. 

To prove uniform boundedness, set g := ||<?||L°°(r); thus the function 
vr := (g) m verifies: 

'Awfl = {x,y)eQ R 

9vr i_ 
-e-^— + P(vr)™ = pg > pg x G Tr 

Vr > (x,y) G Sfl. 

Then, by (|3.3p . ;z# < 5jj in Tr. Thus, the function £r := vr — vr satisfies 
the problem 

{A( r = (x,y)en R 
\Cr>0 (x,y) E^rUTr. 
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The comparison principle then implies Cr > in £Ir, that is 

vr < ||s|ir~(r)> in °fl- 

On the other hand, since z = is a subsolution to problem ()3.5j) . v R > 0; 
(13.41) is then proved. 

To establish in. we consider v R i solution of (|3.1|) in $7^'. Since, by ii., 
v Ri is nonnegative in R > , we have 

'Av R , = {x,y)en R 
dvR ' i / \ — ^ r 

.ur' > (x,y) G 



Then, by (l3~3l) . 



/ (zr - Zi?')+pd.x = 0; 



this implies zr< > zr on T^. To conclude, observe that the function rjR := 
v R f — v R satisfies, 

|A^ = (x,y)eQ R 

\?ir>o (ijjeSflUrjj. 

Hence, by comparison, v R i > v R on fl R . □ 

Next step in the proof of Theorem 13.11 consists in studying the following 
evolutive problem in D,r: 



Aw R = 
w R = 



(x,y)eQ. R , t>0 
(x,y)E'E R , t>0 



dw R 







(w R y< 



(3.6) 



at 



x G Tr, t > 
xeT R , t = 0. 

is a pair of functions (u R ,w R ) 



dy 

(WR = uf 

Definition 3.4. A solution to problem 
such that 

• ur g L°°(T r x (0,oo)) nC^oo);!^^), ur > 0; 

• w R G L°°{n R x (0, oo)), w R > 0, \Vw R \ G L 2 (n R x (0, +oo)); 

• wr\„ , , = vlfi , wr\„ , n . = 0; 

• for anyT > and any i\> G CI'}(?[r~ x [0,T]), V = in dT, R x (0,T) 



o Jo 



(Vwr, V'0} da; dy dt 



pURdtip dx dt 



p(x) ur{x, T)ip(x, 0, T) — uo(x)ip(x, 0, 0) dec. (3.7) 
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Observe that, by standard results, for any g G L°°(Tji) there exists a 
unique solution v G W 1,2 {Q,r) P\ L°°(VIr) of problem 

Av = in Qr 

v = on T,r 

v(x, 0) = g(x) oiiTj}; 

here boundary conditions are meant in the sense of trace. We shall denote 
such a solution by E^g) since it can be regarded as the harmonic extension of 
g to Qr, completed with homogeneous zero Dirichlct boundary condition on 
Sr. Indeed, by standard elliptic regularity results, v G C 2 (f^)nC(f2 fi UE fl ) 
and v(x) = for all x G 



Proposition 3.5. Let assumption JAoJ) be satisfied. Then for any R > 
there exists a solution (ur,Wr) of problem 



Proof. We introduce a time-discretized version of (|3.6[) . Fix any T > 0. 
For any n G N we divide the time interval [0,T] in n subintervals of length 
e = T/n and endpoints [(fc — l)e, fee] for fc = 1, . . . , n. For any k = 1, . . . , n, 



by Lemma 13.31 with g 
the problem 



,fc-l 



'/v 



there exists a unique solution (u R , u>^j ) to 



f a c-fe 

e.fc 
W R = 



dy 



P [ U A 



,k-l 



e,0 



« 



(re, y) G Or 
(x, y) e S fl 

a: G Ttj 
a; G r fl . 



(3.8) 



The solution (u R k ,w R k ) satisfies 



,k-l 



("* ) 



E f 



We can rewrite the problem on Tr in 

6,fe-l 



as 



6,0 



n 



x G r fl , 
x g r fl , 



where ^ : V(A) C L 1 ^) -> £°°(r,R) is the operator defined 



A(v) := - 



1 /SE R (^ 



dy 



with domain 



^):=fe™: ^)el»(r s ), |MU~(r B ) < |MU~(r R )}. 
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The operator A satisfies the following properties. For any e > the mapping 
(I + eA) is one-to-one from T>(A) onto a subspace lZ e (A) C L°°(rn). In fact, 
by Lemma 13.31 for any e > there exists a unique solution of 

Aw e £ k = (x, y) e Q R 

1 0w e ' k 

u 



f_ e L^ = u f-i xeTr (3.9) 



w% k = (x,y)eE R \ 



P dy 

J R ' ' 

in addition, by (|3.3|) . the inverse mapping 

(I + eA)- 1 :K e (A)^L™(r R ) 

is a contraction with respect to the norm of L 1 (Fr). 

A second property satisfied by A is the following rank condition: for any 

e > 

TZ e (A) = L°°{T R ) D L°°(T R ) = V(A) for any e > 0. 

The validity of the two properties above permits to invoke the Crandall- 
Liggett theorem (see [9]) and to infer that u e R converges in L 1 (Fr) to some 
function u R G C([0, T]; L 1 (T R )) , as e -» 0. Furthermore, by ([3^1]) and ( |A^ )- 

(ii)> 

< w^ fe < ||uoIIl°°(r) in n R , uniformly with respect to e, k and i?. 

Then w R k converges weakly-* to some function w R £ L°°(VL R x [0,T]) as 
e 0. 

Using ui^ fc as a test function for p.8j) . integrating by parts, and applying 
Young's inequality we obtain 



< e 



^ W R k 



2 1 

da; dy < 



dx. 



1 + m Jr R 

(3.10) 

Then, adding for k from 1 to n and passing to the limit as e — > we 
discover 



/ / \Vw R fdxdydt< — !— / 



p (u ) m+1 dx. 



Thus, ffljj € L 2 ([0, T]; 7J 1 (Ofl)). Moreover, it is easily seen that w R 
0. Concerning the limit function u R , we observe that, by (|3.10l) . 



E H x(0,oo) 



/) I m r dx < P \ U R I dx < P u o dx, 

for any e > 0; then for e — ^ 0, we get 

/ p(u R (x,t)) m+1 dx< [ pu™ +1 dx, for any t e [0,T]. 
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Now, by choosing appropriate test functions as in |30j and taking into 
account (|3.2[) we have for any k — 1, . . . , n, 



[ ( Vwg fc , V^) dx dy = - I p u R k 1 — u R k tpdx. (3-11) 
Jn R e Jr R 1 J 

Integrating on [tk-i,tk] and then adding over k, we can rewrite the right 
hand side of equality above as 



1 . " . r tk f r 

- 2_. j j P u R^( x >t)~ u R k ( x it) ip(x,0,t) dx dt 



k=i Jt k-i J r B 

7 P uUx,tf ixAt + e) -^ XAt) d X dt 
i Jt r e 

i r f i f T f 

puo(x)tp(x, 0, t) dx dt / / pu R (x,T)'tp(x,0,t) dxdt, 



c 



o Jr R c JT-e Jr 



(3-12) 



where u e R := {w^ 1 , . . . , w^™}. By sending e to 0, (|3.1ip and (|3.12[) give 
(15771 . □ 

In the next Proposition we list some properties satisfied by solutions of 

Bra. 



Proposition 3.6. Lei (ur,wr) be a solution of (|3.6p . TTien £/ie following 
properties are satisfied: 

i. If ur and ur are solutions of (|3.6p corresponding to u$ and uq respec- 
tively, then 

[u R (x,t) - u R (x,t))} + pdx < / [u - u } + pdx; (3.13) 



in particular, if Uq > then u R {x,t) > for every x G Tr and every 
t > 0; 

ii. There exists a constant C independent of R such that < wr < C for 
every x € £Ir and every t > 0; 

Hi. For any R' > R, WR>(x,t) > wn(x,t) for every x £ £Ir and every 
t > 0; 

iv. The inequality 

p(x) [(m — l)tdtWR + mwpt] > 
holds in the sense of distributions in £Ir x (0, oo); 

v. d t u€Ll c ((0,^);L 1 (T)); 
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vi. For any < r < T, 

p(u R ) m+1 (x,T)dx 



[ f \Vw R \ 2 dxdydt+^— f 
Jr Jn R m + 1 Jr R 



1 



m + 1 



p(u R ) m + 1 (x,r)dx; (3.14) 

r 



vii. For any < r < T and any cutoff r/ 6 Cg(f2) with < rj < 1 arte? 



supp(?7) c Bji, 

2 



/ / \Vw R \ 2 rj 2 dxdydt 



m + 1 



p{x) {u R ) m+1 {x,T) V 2 dx 



p(x)(u R ) m+1 (x,T) V 2 dx-2C 2 



771+1 



[ [ \\7n\ 2 dxdydt. 
Jo Jn R 

(3.15) 



Proof. Properties i.—iii. are inherited by {u R , w R ) from those listed in Lemma 
3.31 for solutions to the auxiliary elliptic problem. In particular i. is a conse- 
quence of (|3.3j) . Property ii. follows taking into account that w R is null on 
Ti R and previous property L; property Hi. follows by Lemma 13.31 - Hi. 

Properties iv.-vi. can be obtained arguing as in [101 Sections 3 and 5]. 
Moreover, it is direct to check that 



/ / \Vw R \ 2 n 2 dxdydt+—!— f , 

Jr JQ R m + 1 Jr n 



p{x) (u R y n+i (x,T)r, 2 dx 



1 



m + 1 



p{x)(u R ) m+1 (x,T)ri 2 dx-2 [ r]w R (\7n,\7w R )dxdydt. 
r R Jt Jq r 



Using Young's inequality, we estimate the last term in the computation above 
as follows: 



2 / / i]w R (Vr],Vw R ) dxdydt <2 / / (\Vw R \n)(w R \\Vr]\) dxdydt 
Jo Jn R Jo Jn R 

<\ ( [ \Vw R \ 2 n 2 dxdydt + 2 f f \w R \ 2 \\7n\ 2 dx dy dt. 
2 Jo Jn R Jo Jn R 



The last inequality, together with the fact that, by ii., w 2 R < C 2 , gives 
(I3~15l) . □ 



Proof of Theorem \3.1[ By Proposition 13. 5[ for any R > there exists a so- 
lution (u R ,w R ) of (|3.6p . Moreover, thanks to the uniform boundedness and 
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monotonicity of wr (properties ii. and in. of Proposition 13. 6[) . there exist 
the limits 

lim ur =: u in T x (0, oo), 

R— >oo 

lim Wr =: w in f2 x (0,oo); 

R— ►oo 



again by Proposition 13.61 m.. w e L°°(0 x (0, oo)) and u € L°°(r x (0, oo)). 

We shall prove that (it, w) is a solution to problem (jl.ip . To see this, take 
T > 0, D and V 1 as in Definition 12.1) then select i?o > large enough, such 
that D C Qro- By Definition 13. 4[ for every R > i? we have: 



/ / (Vi«ij, Vi/>) dzdydt + / / pu R d t ipdxdt 
Jo Jd Jo JdDnr R 

p(x) u(x, T)ip(x, 0, T) — uo(x)ip(x, 0, 0) dx. (3.16) 



Furthermore, by Proposition 13.61 w'., since it G L°°(r x (0,oo)) and uq G 
L°°(r), for each compact set K C fi, there exists a positive constant C, 
depending on K but independent of R and T > 0, such that 



/ iVwflJ 2 < C in K. 
o Jk 



By usual compactness arguments, Vu> S L^ oc (i7 x (0, oo)), so w\r is well- 
defined; moreover, letting R — > oo in (|3.16[) we get 



/ ( (Vw, Vip)dxdydt + [ [ pud t i/jdxdt 
Jo Jd Jo Jod 

p(x) u(x, T)ip(x, 0, T) — uo{x)ip(x, 0, 0) da;. (3.17) 



Hence, it is direct to check that (u, w) is also a solution in the sense of 
Definition O □ 

Remark 3.7. By construction, the solution (u, w) given in the proof of The- 
orem [33] is minimal, namely, if (u,w) is a solution to (|2.1I) with u > and 
w > 0, then u > u and w > w. □ 

Remark 3.8. (i) As already pointed out in the Introduction, our construc- 
tion of the solution is slightly different from those in [TU] and [TT]. In fact, 
the strategy in |10j consists in passing first to the limit as R — > oo in the 
discretization (|3.8[) to obtain a solution of an elliptic problem in the whole of 
O, and then in sending the discretization parameter e to zero. Our approach 
consists instead in passing first to the limit as e — > in the discretized elliptic 
problem, in order to get, for any R > 0, a solution of the parabolic problem 
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(13.61) in ttf{. Finally sending R — > oo we get a solution in the whole space. 
A similar strategy has been implemented in [11) , where unbounded cylinders 
Bft x ]R + , instead of bounded domains considered here, have been used 
to invade £1. Note that the use of Sl^ will be expedient in proving Theorem 
14.11 in the following Section. Indeed, we shall construct a suitable family of 
solutions to elliptic problems in bounded domains Sl#. 

(ii) Introducing a varying density p does not bring any additional techni- 
cal difficulty to the proof of existence of solutions to problem (|2.1[) compared 
with those in [TU]- jTT] , in which p is assumed identically 1. 

(Hi) The restriction N = 1 does not play any role yet and all results of 
Section [3] can be easily adapted to show existence of very weak solutions in 
any space dimension. 

(iv) Concerning the regularity required for the initial data, we should re- 
mark that in [10] and [TT] uo is assumed to belong to L l , whereas in [28] and 
[29] - where the (non fractional) porous medium equation with varying den- 
sity has been studied - it is supposed to belong to L^T). As a consequence, 
the solution considered in [10] and [TT] satisfies |Vu>| £ L 2 (Q, x (r, oo)) for 
any r > 0. We assume instead uq nonnegative, continuous and bounded, so 
Mo is only in Lj oc (T) = ^^(Y), but in general uo ^ Lp(T). Hence we can 
infer that \7w £ Lf oc (fl x (0,oo)), as required in our Definition 12.11 but in 
general Vu> <^ L 2 (Q x (t, oo)) . 

(v ) We emphasize that, unlike what have been discussed in (ii) and (iv) 
above, the presence of a variable density and the assumptions on the initial 
data shall entail significant differences in the proof of uniqueness of solutions, 
relative to those in the case in which p = 1 and uo £ L 1 ; see also Remark l4.3l 
Moreover, the hypothesis N = 1 will be essential in our proof of uniqueness; 
see Remark [44] □ 

4 Uniqueness of solutions 

The scope of this Section is to establish the following 

Theorem 4.1. Let assumption QAqD be satisfied. Then problem (|1.1[) admits 
at most one solution (u,w) with u > 0, w > 0. 

To prove Theorem 14.11 we need some preliminary materials. To begin 
with, observe that the function 



e(x,y): 



- log | (x, y) |, (x,y)en 



(4.1) 



satisfies 




(4.2) 
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where do is the delta distribution concentrated at the origin. As a conse- 
quence, (see, e.g., [I] and [5]) for any F € C^(T), F > 0, the function 

w(-,y) := e(-,y)*F 

is, up to additive constants, the unique bounded solution of 

( All = in Q 

{ dU (4.3) 
for a detailed proof of the fact that w solves problem (|4.3p we refer, e.g., to 



The following Lemma will play a crucial role in the proof of Thcorcm l4.1 



Lemma 4.2. Let assumption flApD 6e satisfied, and Rq > be fixed. Let ij) R 
be the solution of the problem: 

[ AU = (x, y) 6 

-— = F x £ r fl (4.4) 
dy 

U = (x,y)eE R 

where R > Rq, F E C°°(T), F > 0, suppf C T Ro . Then the following 
statements hold true: 

(i) For any Ri,R 2 £ (Rq,oo), R\ < R2 we have 

< ip Rl < ipB, 2 in n Rl ; (4.5) 

(ii) There exists a constant M > depending on Rq such that for any 
R > 2Rq we have 

M <^<0 onH R , (4.6) 



i?(logi?-logi? ) - dv R 
v R denoting the outer normal to Yi R . 



Proof. (1) By the maximum principle ?p R > in Q R for any R > Rq; hence 
the function ip Rl — ip R2 (R < Ri < R 2 ) is a subsolution of problem 

'AU = (x,y)en Rl 
dU 

= o x e r Rl 

dy 

U = (x,y)EZ Rl . 
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Then again by the maximum principle we get (|4.5[) . 

(ii) Clearly, by the strong maximum principle it follows that, for any 
R > i?0j 

R < on E fl . (4.7) 



9vr 



Put 



©(•) y)* F =: ipco(x, y), (i,y)en, 

where = 0(|(x, y)|) is defined in (|4.1[) . It is easily seen that for any R > R 
ipoo is a supcrsolution (|4.4|) . By comparison principle, 

< V'oo in fifl. 

Indeed, it could be shown that 

4>oo(x,y) = hm i/j R (x,y), (x,y) £ 0. 
Now, for any R > i?o, consider the problem 



Af7 = 
_dU _ 

dy 
U = 



o x g r fl \r flo 

(x,y) G 
[/ = M (x,y)GS fio , 



(4.8) 



where M := max-^oo. Define 



U R 



Z(x,y) = Z(\(x,y)\) 



M 



e(\(x, y )\)-e(R) 

Q(Ro)-e(R) ■ 



(x,y) £ fl R \ft Ro . 



Using (|4.2[) immediately follows that Z is a supersolution to problem 
for any R > 2Rq. On the other hand, since in particular R > Rq, i/jr is a 
subsolution to the same problem. By comparison we get for any R > 2i?o 

Z>tp R inf2 fl \fi Ro . 

Consequently, since for any i? > 2i? 

Z = i'R = on £r, 

we can infer that, for any (x,y) G £_r, 



dipR dZ _ dZ{R) 



AI 



dv R dv R 



dr 



R(-logR +logR)' 



here r = r{x,y) = \(x, y)\. The above inequality, combined with (|4. 7[) gives 
This completes the proof. □ 
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Now we can prove Theorem l4.ll 



Proof of Theorem \4-l\ Let (u, w) be the solution of problem (jl.lj) constructed 
in Theorem 13.11 such solution is minimal as observed in Remark 13.71 Let 
(w, w) be any other solution of the same problem and set 

K := ||H|i°o ( nx(o,oo)) ■ 

We claim that 

/ / [u m Fdxdt = (4.9) 

Jo Jr 

for any T > 0, F G C °°(F). 

By the Claim the conclusion follows. In fact, in view of the arbitrariness 
of F, (|4.9p implies 



u = u in F x (0,T), (4.10) 

whence the conclusion. 

It remains to prove the Claim. To this aim, without loss of generality, we 
suppose supp(F) C r_R for some R > 0, F > 0, F ^ 0. Take R > 2R and 
let ipn(x, y) be the solution of problem (|4.4j) . Using tpR as a test function for 
the equation (|2.1[) (see Definition I2.1[) we obtain, for any r > 0, 



w db dt 
o Jdn R ovr 



p(x) u(x,t) — uo(x) ipu(x,0)dx 



^A rU 

w—-dx dt 
o Jr R oy 



w dx dy dt 
o J?. R ovr 



Jt 



p(x) u(x,t)—uq(x) iPr(x, 0)dx 

R 

u m Fdxdt - [ [ w-^dxdydt 
Jo Js n our 

p(x) u{x,t)—uq{x) iPr(x,0) dx. 

R 

Similarly, using ipR as a test function for u, we get 



(4.11) 



lo Jr 



u m Fdxdt 



^ A A Af 

w—^-dx dy dt 

o Js R ovr 

p(x) u(x,t)—uo(x) iPr(x, 0)dx. (4.12) 
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o Jt 



[u m -u m ]F(x) dxdt 



By subtracting (|4.12l) to (|4.11|) we have 

p(x) [u(x, t) - u(x, t)] ip R (x, 0) dx 

= - I I {w-w}^dxdt. (4.13) 
Jo Jy, r dv R 

Since F > 0, i/jr > and, by minimality of u, u > u, equality (|4.13[) with 
t = T gives 



T 

Jt, 



l ]F(x) dxdt 

[ [ { w — w\ dx dt 

Jo Js R OVr 

By recalling that |Er| = irR/2 and by using (|4.6|) we obtain 



< liminf 

R— >-oo 



(4.14) 



/ / {iy-w}^^da 
Jz R OVR 



dt 



< KttRsvlp 

< K-k- 



dipR 



Qvr 
M 



0, as R — > oo. 



(-logiJo + logi?) 

Hence, the right hand side of (|4.14[) is zero and the claim (|4.9[) follows. This 
concludes the proof. □ 

Remark 4.3. Uniqueness results in [10] and [TT] cannot be applied to prob- 
lem (jl.lj) . as they require p = 1. Furthermore, even in the case in which 
p is identically 1, the results in [TU] and [TT] cannot be adapted to prove 
uniqueness of solutions in the sense of Definition 12.11 □ 



Remark 4.4. If N > 2, a result similar to Lemma [4.21 could be proved. In 
this case, the estimate (|4.6p must be replaced by the following: 



M 
R^ 



<^<0, 

OVR 



on Y,r, 



(4.15) 



for some M > 0, for any R > 2Rq. Hence we cannot get the conclusion in 
the previous proof. A similar issue arises for problem (|1.3j) when N > 3; 
however, methods used in [24] to get the conclusion in that case cannot be 
adapted to the present nonlocal situation. □ 



Remark 4.5. Note that the proof of uniqueness for problem (|1.3|) when 
N = 1,2 is quite different from the previous one (see [T7], [27]). However, 
the estimate (j4.6[) is crucial in [TT] for N = 2, too. Observe finally that the 
hypothesis p 6 L oa (M 2 ) used in pTJ is not required here. Thus our arguments 
can be adapted in order to prove uniqueness for problem (|1.3j) in the case 
N = 2 and p possibly unbounded. □ 
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